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Consumption of traditional aniseed alcoholic beverage, raki, adulterated with methanol results in deaths, 
therefore, its detection is an important issue. In this study, mid-infrared spectra of pure and methanol 
adulterated (0.5-10Z (vol/vol)) raki samples were collected with an attenuated total reflectance attach- 
ment of a Fourier-transform infrared spectrometer. Principal component analysis was used to discrimi- 
nate pure and adulterated raki samples, then, a partial least sguare model was constructed to 
determine the adulterant methanol content in raki using mid-IR spectral data. A minimum threshold 
level of 0.55 methanol in raki samples was successfully detected. A good prediction model for determi- 
nation of methanol adulteration ratio in raki samples was also constructed (R2 - 0.98 and RPD - 8.35). 


© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 


Aniseed flavoured distilled spirits are produced in different 
Mediterranean countries and are known with various names such 
as anis, arak, pastis, ouzo, sambuco and zivania (Anli & Bayram, 
2010). Raki, a traditional Turkish alcoholic beverage, also belongs 
to this group. It is the second most produced alcoholic beverage 
in Turkey (20,083 L absolute alcohol in 2009 (TAPDK, 2010)) and 
is described in Turkish Food Codex (2005) as a spirit that must 
be produced in Turkey by double distillation of suma or suma 
mixed with agricultural-based ethanol and flavoured with aniseed 
(Pimpinella anisum). Suma is produced mainly from raisins, molas- 
ses and/or grape must and is a distillate with a maximum 94.52 
ethanol content. Suma can be also mixed with ethanol from an 
agricultural origin and this final mixture must contain at least 
65X suma. Suma or suma-ethanol mixture is distilled in copper 
alembics with a maximum capacity of 5000 L and three groups of 
distillates, head, medium and tail, are separated (Anli & Bayram, 
2010). Medium distillate is used for raki manufacture and diluted 
to 45 alcohol. After the addition of sugar product is stored for at 
least one month before filling. 

Alcoholic beverages could be mixed with ethanol, methanol, 
water and aromatic compounds to gain illegal profit (Pontes 
et al., 2006). Consumption of the excess amounts of some of these 
adulterants could cause severe health problems. Severe metabolic 
disturbances, blindness, permanent neurologic dysfunction and 
death are the results of methanol ingestion. Although methanol it- 
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self is harmless, its metabolites are toxic and its consumption 
causes methanol poisoning in both humans and animals. Methanol 
spreads rapidiy into the body via all routes. After dispersion in 
body compartments, methanol is firstliy digested into formalde- 
hyde by alcohol dehydrogenase enzyme. Formaldehyde is thirty- 
three times more toxic than methanol, and it is speedily converted 
into formic acid by aldehyde dehydrogenase enzyme. Finally, for- 
mic acid is metabolised into carbon dioxide then water is moved 
away from body (Kapur, Vandenbroucke, Adamchik, Lehotay, & 
Carlen, 2007). Lethal dose of methanol in humans is related to 
age, sex, metabolic rate and methanol blood level. It was estimated 
that consumption of 50-75 g of methanol causes deaths due to 
methanol intoxication; however, there are also reports showing 
that toxication level could be aslow as 11.5 g (Di Maio & Di Maio, 
2001; Fidan, Denli, & Anli, 1996). Many death cases related to 
methanol intoxication due to consumption of fake raki have been 
reported in Turkey. According to literature, the total number of 
deaths due to methanol poisoning which was arisen from con- 
sumption of adulterated raki, has been more than 300 from 1990 
to 2011 (Eke, Buyuk, Dinc, & Citici, 2007; Gulmen, Meral, Hilal, 
Akcan, & Cekin, 2006; Yaycı & Inanici, 2005). Generally, chromato- 
graphic methods are applied for the detection of adulteration of 
alcoholic beverages (Savchuk et al., 2001; Yilmaztekin & Cabaroglu, 
2011). However, there has been an increasing interest towards the 
development of rapid technigues such as the use of different spec- 
troscopic methods. Front face fluorescence spectroscopy was used 
in differentiation of brandies and wine distillates (Sadecka, 
Tothova, & Majek, 2009). Zivania, a traditional Cypriot aniseed fla- 
voured distilled beverage, was characterised by NMR spectroscopy 
and the data was used satisfactorily to differentiate zivania from 
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other similar beverages (Petrakis et al., 2005). Distilled beverages, 
whiskey, brandy, rum and vodka were classified successfully using 
near infrared spectroscopy and verification of the mixing of these 
beverages with 5-105 of ethanol, methanol and water was also 
achieved with the same technigue in combination with chemomet- 
ric analysis (Pontes et al., 2006). 

Mid-infrared (mid-IR) spectroscopy in combination with che- 
mometric technigues has become a very effective tool in detecting 
the adulteration of various food products such as honey, different 
types of oils, coffee and fruit juices (Briandet, Kemsley, & Wilson, 
1996; Gurdeniz & Ozen, 2009; Kelly & Downey, 2005; Sivakesa & 
Irudayaraj, 2001). There are also examples of the application of 
mid-IR spectroscopy in literature for authentication of wines from 
different origins as well as detection of its mixing with adulterants 
such as giycerol (Dixit, Tewari, Cho, & Irudayaraj, 2005; Tarantilis, 
Troianou, Pappas, Kotseridis, & Polissiou, 2008). 

This study aims at detecting and guantifying the adulteration of 
raki with methanol with mid-IR spectroscopy in combination with 
chemometric technigues. 


2. Materials and methods 
2.1. Samples 


Raki samples were obtained from local markets and cover most 
of the brand names and their sub-brands produced in Turkey. More 
than one sample from some brands was also included into the 
sample set. Two ouzos, aniseed drink from Greece, were also tested 
as samples. Totally, 34 samples were analysed. Samples, their raw 
materials and process types according to descriptions on the labels 
and information obtained from producers are listed in Table 1. The 
raki samples were stored in their original bottles at 4 *C. The lids of 
the bottles were covered further with parafilm in order to avoid 
alcohol volatilisation. 

Six representative raki samples (1, 10, 15, 16, 17, 24) were se- 
lected and they were adulterated with 0.5-10Z (v/v) methanol 
(99.95, Sigma). Spectra of all pure and adulterated samples were 
obtained with a Fourier-transform infrared (FT-IR) spectrometer. 


Table 1 
List of raki samples. 


Sample number Oo Production process Oo Grapetype Oo Alcohol content (X) 


1 Distillated twice Fresh-dried o 45 
2 Distillated twice Fresh-dried o 45 
3 Distillated twice Dried 40 
4 Distillated twice Fresh 45 
5 Distillated twice Dried 45 
6 Distillated twice Dried 43 
7 Distillated twice Dried 45 
8 Distillated twice Fresh-dried o 45 
9 Distillated twice Dried 

10 Distillated twice Dried 45 
11 Distillated twice Dried 45 
12 Distillated twice Dried 45 
13 Distillated twice Fresh-dried o 43 
14 Distillated twice Fresh 50 
15 Distillated thrice Dried 47 
16 Distillated twice Dried 45 
17 Distillated twice Fresh 47 
18 Distillated twice Dried 45 
19 Distillated twice Dried 

20 Distillated thrice Dried 45 
21 Distillated thrice Fresh 45 
22 Distillated twice Dried 43 
23 Distillated twice Dried 43 
24 Distillated twice Dried 45 
25 Distillated twice Dried 45 
26 Distillated twice Dried 45 
27 Traditional ouzo 38 
28 Traditional ouzo 38 


2.2. Mid-IR analysis 


All infrared spectra (4000-650 cm!) were obtained using a 
Perkin Elmer Spectrum 100 FT-IR spectrometer (Perkin Elmer 
Inc., Wellesley, MA) with a deuterated tri-glycine sulphate (DTGS) 
detector. Samples were placed on a horizontal attenuated total 
reflectance (HATR) sampling accessory with a ZnSe cıystal 
(45 deg. Trough Plate) and covered with a lid to prevent evapora- 
tion of alcohol. The number of scans for each spectrum was 64. 
The resolution of the instrument was set to 4cm |, ZnSe crystal 
was cleaned with ethanol and water in between sample runs. At 
least three measurements were collected for each sample. 


2.3. Multivariate analysis 


Spectral data obtained with an FT-IR spectrometer were ana- 
Iysed by using multivariate statistical technigues with SIMCA soft- 
ware (SIMCA P-10.5 Umetrics Inc. Sweden). Whole spectrum 
(4000-650 cm!) and also three useful spectral intervals selected 
to avoid interferences and noises were used in the data analysis. 
These selected spectral regions were 1565-965, 1900-1700 and 
3040-2800 cm”! (Lachenmeir, 2007; Versari, Parpinello, Scazzino 
& Rio, 2010). 


2.3.1. Pre-treatment of data 

Spectral data set was transformed into suitable form by pre- 
processing technigues of scaling and centering. Raki samples con- 
stitute the columns and the variables (spectral data) are located in 
the rows. Pre-processing can be improved by using mathematical 
procedures including signal corrections, some various derivative 
forms and using some filtering technigues such as wavelet com- 
pression of spectra (WCS), wavelet in combination with orthogonal 
signal correction (W0OSC) and orthogonal signal correction in com- 
bination with wavelet (OSCW) (Casela, Sinelli, Oliveri, Valentina Di 
Egidio, & Silvia Lanteri, 2010). For obtaining better classification 
with principal component analysis (PCA), wavelet compression of 
spectra (WCS) was preferred as spectral filtering technigue in this 
study. Daubechies-10 function wavelet was chosen and confidence 
interval was set to 99.95. De-noise and sophisticated signals are 
compressed with WCS technigue which depends on retaining the 
significant coefficients from the representation of the data in the 
new ordinate system (Eriksson, Johanson, & Wold, 2001). For par- 
tial least sguare (PLS) analysis, OSCW technigue was applied. 


2.3.2. Classification 

Spectral data of pure and methanol adulterated raki samples 
were analysed by PCA for classification. The data matrix consists 
of observations represented by raki samples and variables repre- 
sented by FT-IR spectral data. The classification method is used 
to determine the group of raki samples based on the similarities 
among members of the same data class (Pontes et al., 2006). 
Assignment of the sufficient number of principal components 
(PCs) depends on the goodness of fit, which is represented by the 
parameter R?, and predictive ability of the model (Eriksson, Johan- 
son, & Wold, 2001). 

Score and Coomans' plots were constructed as part of classifica- 
tion analysis. When a grouping was observed in score plot then 
Coomans' plot was also used. The Coomans' plot compares the dis- 
tance to the model results against the distance from the model 
centre for unknown samples of selected models (Esbensen, Guyot, 
Westad & Houmoller, 2002). Coomans' plot comprises of four re- 
gions, two regions specify pure and adulterated classes defined 
by the models, common region indicates overlap of the two classes, 
and the other region indicates the observations far from two 
classes. 
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2.3.3. Ouantification 

PLS was employed to determine adulteration rate of methanol 
adulterated raki using both whole and selected spectral ranges. 
PLS analysis was performed to establish a regression model for pre- 
diction of adulterant concentration (variable matrix Y) using FT-IR 
spectra (variable matrix X). 

Obtained data set was randomly separated into two groups as 
calibration (2/3 of samples) and validation (1/3 of samples) sets. 
The predictive ability of the model was expressed by some param- 
eters and visualised with prediction plots of constructed models. 
These parameters are root mean sguare error of calibration 
(RMSEC), root mean sguare error of prediction (RMSEP), relative 
error of prediction (REP), residual predictive deviation (RPD) and 
the regression correlation coefficient (R2) both for calibration and 
validation models. RMSEC and RMSEP values are used to evaluate 
the performance in the prediction process. RMSEC refers to the cal- 
ibration uncertainty that can be expected for predictions. A good 
model would have small value of RMSEC and RMSEP. Assessing 
the accuracy of the model can be also expressed using RPD value. 
Generally evaluating all these parameters gives an idea about the 
predictive efficiency of the model (Esbensen et al., 2002; Tamaki 
& Mazza, 2011). RMSEC, RMSEP and RPD are defined with follow- 
ing eguations: 


RE ARENEN 
27 — yp)” 
RMSEC— Yİ. 
n—2 
n 
Yy) 
RMSEP — YE! 


where n is the number of samples used in each set; SD is standard 
deviation in each set, y; is the actual value for i-th sample, f; is the 
predicted value determined by FT-IR for the same sample and y is 


2.01 


the mean of each set. Criteria for R?, RPD and slope used to evaluate 
the performance of predictive models are listed by Tamaki and 
Mazza (2011). 


3. Results and discussions 


Fig. 1 shows the typical spectra of the raki samples obtained in 
this study. All pure samples have similar spectra; however, there 
are changes in the intensity of some peaks depending on the sam- 
ple. These differences are mainly for peaks at 3300, 2975-2900, 
1640 and 1450-1250 cm |, FT-IR spectra show absorption bands 
at different freguencies and these bands are attributed to the var- 
ious functional groups. The two distinct absorption bands around 
3600-3200 and 1700-1565 cm” are associated with O-H stretch- 
ing vibrations. The peak at 3040-2800 cm”! is due to the stretch- 
ing vibratton of C-H bond. The peaks between 1565 and 
965 cm! are in the fingerprint area where many different IR 
bands, including those corresponding to the vibrations of the C- 
O, C-C, C-OH, C-H, and CEN bonds occur. Bands at 2960- 
2850 cm! (C-H stretch) and 1083-1045 (C-O stretch) cm! are 
attributed to ethanol while various groups belonging to organic 
acids, carbohydrates and phenols provide peaks at fingerprint re- 
gion of 1500 cm! and below. SpecificalIy, it was reported that or- 
ganic acids have peaks at 1150-1060cm ! (C-O stretch) and 
1420-1320 cm! (0-H bend) and peaks around 1100-900 cm! 
belong to residual carbohydrates (Polshin et al., 2011). The region 
from 4000 to 3600cm ! and the spectral range from 965 to 
650cm! contain very little useful information (Fernandez & 
Agosin, 2007; Lachenmeier, 2007). Therefore, in addition to whole 
spectral range, a selected range (1565-965, 1900-1700 and 3040— 
2800 cm) was also used in statistical analysis to reduce noises by 
eliminating areas with no peak and excluding -OH peaks which 
were very dominant. 

The firstaim of this study is discrimination of methanol adulter- 
ated and pure raki using PCA. Adulteration range that was used in 
this study was 0.5-10Z (vol/vol). Raki could be produced from 
fresh, dried or fresh and dried grapes. Either suma alone or suma 
and ethanol from an agricultural source could be used in distilla- 
tion. Also, distillation could be applied twice or three times during 
production. Six raki samples which are representative of different 
raw materials and processing technigues were adulterated so that 
it would be possible to take into account of the differences be- 
tween raki samples. IR spectra of methanol adulterated and pure 
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Fig. 1. Typical raki (gray line) and methanol adulterated (black line) mid-infrared spectra of raki. 
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raki are provided in Fig. 1. As can be seen from spectra, distinctness 
arising from methanol content could be determined by comparing 
several regions. The first big difference corresponds to region 
around 3500 to 3200 cm ! that can be attributed to O-H stretch- 
ing. Another difference in region 2900-2800 cm! arises from 
C-H stretching and its source could be methyl groups. The other 
distinction between pure and adultered raki is in the fingerprint 
area (between 1565 and 965 cm |). The differences between pure 
and adulterated sample spectra could be attributed to the absor- 
bance bands specifically to methanol (spectrum not shown). 

Both the whole and selected ranges of spectra were used for 
classification. WCS filtering was also applied to both cases. Accord- 
ing to score plots (not shown), similar results were obtained using 
both the whole spectra and the selected range. Adulterated raki 
samples generally located on the left side of score plot while pure 
ones at the right side. There are some adulterated samples which 
are located in pure samples. When selected range of spectra was 
the data set 8 PCs were used to construct the model and 2 PCs ex- 
plained 91.105 of the variance. 

Since a separation was observed between pure and adulterated 
samples further analysis of data was performed with SIMCA class 
models and Cooman's plot is constructed. Cooman's plot for the 
classification of pure and adulterated raki samples shows good 
separation of adulterated raki with the ratio of 0.54 to 10 from 
pure samples (Fig. 2). By using whole spectral range (4000- 
650 cm !), 12 PCs and 16 PCs were used to construct class models. 
Pure raki samples are located in the upper side of the plot and sep- 
arated from adulterated raki samples and there were no samples in 
the common region. In addition, samples also show grouping be- 
tween each other according to their adulteration ratio. Detection 
of adulteration in alcoholic beverages was investigated in other 
studies. In the study of Pontes et al. (2006); 69 samples of alcoholic 
beverages (whiskey, brandy, rum and vodka) were investigated to 
detect adulteration with water, ethanol or methanol by using 
near-infrared spectroscopy and chemometric methods. Successful 
results were also obtained for the verification of adulteration as a 
result of this study. 

FT-IR in association with multivariate chemometric technigues 
was also employed to guantify the amount of methanol in adulter- 
ated raki samples. Results indicated that PLS model based on IR 
data was proven suitable as a practical analytical method for pre- 
dicting adulterant content in raki samples in the volume fraction 
ranged from 0.5 to 10. Fig. 3 shows observed vs. predicted values 
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of methanol concentration in raki for calibration and validation 
sets, RPD, slope and R? values were calculated for calibration and 
prediction sets to test the accuracy of the model. The accurate 
and similar prediction results were obtained by PLS using both 
whole and selected spectral ranges; for whole spectral range, 
Reaı2 > 0.9838, o Rya2 x0.9914, ORMSEC-0.44, RMSEP-0.35, 
RPD > 8.35 and slope - 0.9 values were obtained and for selected 
range Rcaı2 > 0.9865, Ryaı2 > 0.9874, RMSEC 0.40, RMSEP - 0.38, 
RPD > 8.29, slope > 0.99 values were calculated. According to crite- 
ria for RPD and slope values given by Tamaki and Mazza (2011), R? 
values greater than 0.9 and RPD values greater than 3 provide 
excellent predictions and slope between 0.9 and 1.1 means a very 
reliable model. Therefore, PLS regression of FT-IR spectral data pro- 
vides excellent prediction and a very reliable model for the deter- 
mination of methanol adulteration ratio in raki samples. Besides, 
small RMSEC and RMSEP values and small differences between 
these values indicated that the created models provide satisfactory 
prediction. 

Mid-IR spectroscopy in combination with chemometrics tech- 
nigues is successful in detection and guantification of methanol 
adulteration of a Turkish aniseed alcoholic beverage, raki and could 
be used as a rapid analysis tool especialIy for control laboratories of 
regulatory agencies. Lower detection limit of adulteration with this 
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Fig. 3. OSCW PLS regression of actual vs. predicted adulteration of calibration and 
validation sets (whole spectral range: 4000-650 cm |). 
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Fig. 2. Cooman's plot for the classification of pure and adulterated raki samples using spectral FT-IR data (pure samples and adulterated raki samples are represented with 


black and gray colours, respectively) (whole spectral range: 4000-650 cm |). 
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method is 0.52 (vol/vol) while tolerance of people to methanol 
intoxication could be aslow as 11.5g. 
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